Effect of charge on the dynamics of an acoustically forced bubble 
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The effect of charge on the dynamics of a gas bubble undergoing forced oscillations in a liquid 
due to incidence of an ultrasonic wave is theoretically investigated. The limiting values of the 
possible charge a bubble may physically carry are obtained. The presence of charge influences 
the regime in which the bubble's radial oscillations fall. The extremal compressive and expansive 
dimensions of the bubble are also studied as a function of the amplitude of the driving pressure. 
It is shown that the limiting value of the bubble charge is dictated both by the minimal value 
reachable of the bubble radius as well as the amplitude of the driving ultrasound pressure wave. 
A non-dimensional ratio C, is defined that is a comparative measure of the extremal values the 
bubble can expand or contract to and find the existence of an unstable regime for C, as a function 
of the driving pressure amplitude, P s . This unstable regime is gradually suppressed with increasing 
bubble size. The Blake and the upper transient pressure thresholds for the system are then discussed. 

PACS numbers: 43.35.Ei, 43.25.Yw, 43.35.H1 
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I. INTRODUCTION 



The study of bubble dynamics and cavitation has a 
long and interesting history in the scientific literature. 
One of the earliest works was that of Lord Rayleigh PJ 
in his study of cavitation phenomena, motivated by 
the need to understand and minimize the damage to 
ships' propellors due to cavitation (the low pressure 
on the surface of the propeller blades causes the liquid 
in contact with the surface to spontaneously form 
unstable bubble clouds which often self-organize into 
dendritic structures [2(- These bubble clouds implode 
with enormous force resulting in serious damages to the 
propeller). 

Cavitation, bubble formation and dynamics are present 
in different instances and situations. Analyses and stud- 
ies of the phenomena have been motivated by and have 
explained very distinct natural, practical phenomena. 
Small amplitude oscillations of a gas bubble in a liquid 
were studied by Minnaert [3|- His work showed the 
important contribution of radial oscillations of entrailcd 
air bubbles in the sound heard from running water. 
In nature, the snapping shrimp uses rapid closure of 
its claws to generate cavitating bubbles which stun its 
prey [J|. Bubble formation and kinetics contribute to 
fluid flow in biological systems in blood and cells in the 
micron scale [5| . In the presence of incident ultrasound 
waves, bubbles can also enhance the rate of chemical 
reactions [y, [7[ . Bubble formation and cavitation can be 
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recreated under controlled conditions in the laboratory 
by subjecting a liquid in a container to a standing 
ultrasonic wave, setting up a pressure field within the 
liquid. When the driving pressure amplitude of the 
sonic field becomes larger than the ambient pressure, 
the pressure in the liquid becomes negative. When this 
pressure exceeds the vapour pressure of the liquid, local 
evaporation is caused. The liquid 'breaks' up forming 
tiny micron size cavitation bubble clouds, which implode 
violently within a very short period of time. Frenzel 
and Schultes demonstrated that these cavitation clouds 
emitted low intensity visible light j8|. This phenomenon, 
wherein light is emitted by a gas bubble in a liquid due 
to its rapid expansion and collapse when ultrasound 
is incident on it, is known as sonoluminescence. This 
has also contributed to an extensive study of bubble 
dynamics in the context of sonoluminescence. Several 
other studies have followed in the literature, including 
those of Gaitan and others [9H12J . The radial oscillations 
of gas bubbles in a liquid that are caused due to incident 
acoustic waves cannot be described trivially. These 
are a type of driven nonlinear oscillations that can 
greatly depend on initial conditions and can be chaotic 
in nature. This has been shown conclusively (see, for 
example, [13l4l7| and references therein). Reviews may 
be found in, for example, [3, Il3l.[l8j. 



The presence of electric charge on bubbles in fluids has 
been reported in the experimental literature. [19l423l |. 
In this paper we use the Rayleigh-Plesset equation mod- 
ified by Parlitz, et al., further modified to take into ac- 
count the presence of charge on the bubble, and study 
the effects of charge, driving frequency and amplitude of 
the ultrasonic driving field on the dynamics of the bub- 
ble, with the specific heat ratio being taken as 5/3 in our 
calculations, consistent with the gas in the bubble being 
a monatomic ideal, inert gas under adiabatic conditions. 



In Section II, we describe the model used; the natural 
frequency of oscillation of the bubble is obtained and the 
phase plots show that the maximal radial velocity for 
the charged bubble greatly exceeds that of the neutral 
bubble. In Section III the Blake radius and threshold 
for the charged bubble are calculated. This is followed 
by a discussion of Raylcigh collapse and the influence of 
driving pressure amplitude in Section IV. We show that 
for a bubble of given ambient radius Rq, there exists a 
cut-off value for the maximum charge it can carry, which 
is dictated by the van der Waals hard-core radius. In 
Section V we introduce a new dimensionless ratio ( as a 
measure of the maximum and minimum radius a bubble 
can achieve under acoustic driving. This when plotted as 
a function of pressure amplitude clearly captures the po- 
sitions of the Blake threshold and the upper critical tran- 
sient pressure threshold for acoustic cavitation, and their 
distinct dependence on the pressure amplitude, driving 
frequency and bubble charge. Bubbles that are present 
in various systems can carry a non-zero charge. Hence, a 
proper understanding of their behavior and dependence 
on pressure, driving frequency, ambient radius and other 
parameters, is essential. This is the subject of investi- 
gation of this paper. The results obtained in our work 
therefore become very useful in the context of medical 
uses of acoustic cavitation and ultrasound for diagnostic 
and therapeutic purposes. 
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FIG. 1: Plots of a bubble's radius (above) & radial velocity 
(below) as a function of time (P s = 1.35Po, Ro = 2/j,m, v = 25 
kHz). The curve for the charged bubble (with charge 0.1 pC) 
reaches a larger magnitude (for both radius and velocity plots) 
than for the neutral bubble. (Color Online). 



II. THE MODEL 

We shall briefly describe the model and the set of equa- 
tions being used to describe the radial oscillations of gas 
bubbles in liquids. Various models have investigated dis- 
tinct, disparate limits of bubble-collapse. The original 
work of Rayleigh assumed the surrounding liquid to be 
inviscid and incompressibile [l[ . Plesset and others [2J- 
[26j have included viscosity, surface tension etc. Keller 
and Kolodner used the same expression but with a mod- 
ification for accounting for acoustic radiation by the bub- 
ble by considering the liquid as slightly compressible [27| ■ 
Keller and Miksis further included all these modifications 
- of viscosity, surface tension, the incident sound wave 
and acoustic radiation, in one model to obtain a modi- 
fied equation [281 ] . 

The model we use to describe the system is based on 
the earlier equations of Prosperetti, Parlitz and Keller 
and Miksis and others [T7|, I28I430J and which are modifi- 
cations of the Rayleigh-Plesset equations [l|, [24|, [25[ ■ 

The system wc consider is that of a bubble suspended 
in a liquid and that has some net charge. The presence 
of charge on a bubble is not speculative. As mentioned 
earlier, electric charge has been found on bubbles under 
acoustic forcing. Various reasons could cause charge to 
be present on the bubble, e.g., the migration of ionic 
charges in the liquid onto the bubble surface, although 
the exact mechanism has been debated. See, for exam- 
ple, the work of Alty [H HJ and Akulichev Q. That 



gas bubbles in water can carry a charge has been clearly 
demonstrated long back [20, [2l[. Indeed, it had been 
shown earlier [lj| that air and other gas bubbles in water 
become negatively charged. More recently, this has been 
shown by Shiran and Watmough (23[. Therein, bubbles 
placed in an electric field are clearly demonstrated to 
veer towards one of the electrodes. Other work on the 
presence of charge on bubbles in different situations 
include that done by Bunkin and Bunkin [31| . 
The study of the dynamics of charged gas bubbles in 
fluids is important because such systems find various 
applications, including in medicine - a widely used 
one being in the medical application of ultrasound. 
While it has been overlooked in several models, various 
investigations of charged bubbles have appeared in 
the literature in different contexts. These include, for 
example, [32j. Their work is however of very limited 
scope since they consider the specific case for which 
the polytropic index T takes the value 4/3 and for this 
case it turns out that the terms containing the charge 
mutually cancel out in the Rayleigh-Plesset equation. 

We take V = 5/3 consistent with taking the heat trans- 
fer across the bubble to be an adiabatic process. For the 
sake of simplicity, we assume the charge to be strictly 
limited to, and uniformly distributed on, the surface of 
the bubble. 

The surface charge density of the bubble thus changes as 
the bubble expands and contracts along with the acous- 



tic forcing being externally applied to it. As mentioned 
above, we use a modified form of the Rayleigh-Plesset 
equation, one introduced by Parlitz, et al. [13], which is 
equivalent in first order in (1/c) to that of Keller and 
Miksis [28[ (c being the velocity of sound) and which in- 
cludes, approximately, the sound radiation which is the 
most significant contribution to damping at higher am- 
plitudes of bubble excitation. In the presence of charge, 
this equation describing the dynamical evolution of the 
bubble radius R in time gets further modified. 

Considering the charged bubble as a non-conducting 
charged shell with a constant charge Q, the electrostatic 
pressure can be calculated as has been done by Akulichcv, 
Atchley and others [22|, l3ll - l33f . The presence of charge Q 
(or surface charge density S) causes the inclusion of an 



electrostatic pressure term 2nS 
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where Rq is the ambient bubble radius, Po is the static 
pressure of the, P s and lo = 2-kv denote respectively 
the amplitude and angular frequency of the driving sonic 
field, P v = 2.34kPa is the vapour pressure of the gas, a, p 
and r\ denote respectively the surface tension, density and 
viscosity of the liquid surrounding the bubble, c is the 
velocity of sound in the liquid. In this work, for the pur- 
pose of the numerical results reported, we consider water 
to be the liquid, with p = 998kg/m 3 , 77 = 10~ 3 Ns/m' 2 , 
c = 1500to/s, Po = WlkPa, a = 0.0725N/m, e = 85e 
where eo is the permittivity of vacuum. 
The presence of charge Q modifies the influence of sur- 
face tension, reducing its effective value (the effective sur- 
face tension changes from a for the uncharged case to 
a— Q 2 /(16ireRo)) and induces several interesting changes 
to the dynamics of bubble oscillations. 
When a pressure wave is incident on a bubble in a liquid, 
the difference in pressure can cause expansion and rapid 
collapse of the bubble, this being followed immediately 
after by further, smaller oscillations which are termed as 
afterbounces. This entire sequence of a maximal expan- 
sion of the bubble followed by collapse and afterbounces, 
is repeated in each cycle when we have sinusoidal forcing 
by ultrasound. When a bubble is in the sonoluminescent 
regime, light emission occurs shortly after the bubble's 
violent contraction to a minimal radius, before the onset 
of aftebounces and repetition of this sequence. Details 



of the phenomenon of sonoluminescence may be found 
in, for example, the review articles [lj,[3a]. We recast 
eqn.(l) in dimensionless form for ease of evaluation by 
redefining the radius R through r = R/Rq and time t 
through r = tut, r and r being the new dimensionless ra- 
dius and time variables. Using the static pressure Po as 
the reference pressure we also define the dimensionless 
quantities , P* v = P v /Pq and P* s = P s /Pq. Using an 
ovcrdot to now denote differentiation with respect to r, 
the dimensionless time variable (rather than with respect 
to £), eqn.(l) can be written in dimensionless form as 

r \ f 2 ( f 

1 )rr + Fr + —\3 

c»/ 2 \ c* 

G (1 - P, v + M) (±\ fl + ^(l-3T) 

r 4 \ c* J r \r 
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are all dimensionless constants. 

This form is used when solving the equation numerically. 
In what follows, we will use the dimensional form of the 
equation everywhere. It should be understood though, 
that data points shown in all the graphs have been 
obtained after numerical evaluation of the corresponding 
dimensionless quantities, followed by rescaling by the 
appropriate multiplicative factors to obtain the variables 
in physically realizable units. 

Figure [TJ shows a plot of the radius of the bubble as a 
function of time, obtained by solving equation ([2]), for 
charge present on, as well as absent from, the bubble 
surface. 

For a given driving frequency of the pressure wave, the 
maximum radius attainable by the bubble increases 
with charge present. Conversely, the minimum radius 
achievable also reduces with increasing charge. These 
are as expected. 

The difference in behaviour can be seen more emphati- 
cally in a plot of the bubble's radial velocity as a function 
of time (in Figure [TJ) . In the plot shown, in the presence 
of charge, the maximum velocity increases to more than 
five quarters of its uncharged value. This gives us a 
picture of the dynamics consistent with the time-series 
of the bubble radius, namely, that the bubble oscillations 
become more violent in the presence of charge, causing 



the bubble to expand to a greater extent and contract 
to a lesser minimal radius than in the absence of any 
surface charge. 

The phase portrait for the system is shown in FigJ2]for 
both a charged as well as uncharged bubble. The phase 




FIG. 2: R vs. R phase plots for charged and uncharged bub- 
ble. The charged case (Q — 0.1 pC) spans larger magnitudes 
of velocity and radius. Ro = 2/im, P s — 1.35Po,f = 25kHz. 
(Color Online). 



linearising it, we obtain: 



x + j3x + uj'qX = 0, 



(4) 



where the damping coefficient /3 and the natural fre- 
quency of the oscillator are given by 
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where only terms linear in x and its derivatives have been 
retained. In eqns.(5), the quantity <j) defined as 
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is the equilibrium gas pressure in the bubble. 
The solutions to the damped equation eqn.(4) represent 
small amplitude oscillations. Substituting the values of 
the various parameters in eqns.(5) yields a natural fre- 
quency ujq for a micron-sized bubble that is in the MHz 
range, in conformity with observations in the literature. 



plot for the charged bubble is larger than that of the 
uncharged case, reaching larger values of the magnitudes 
of both radius and velocity. One can intuitively expect 
this. During the expansion part of the cycle, the charges 
present on the bubble surface move away from each 
other as the bubble expands, decreasing the surface 
charge density of the bubble wall. The change in the 
value of the maximal radius of the bubble due to the 
presence of charge, though present, is small. Due to the 
decrease in surface tension, we would expect the radius 
of the bubble to grow slightly more than in the case 
when there is no charge present. On the other hand, 
during the collapse, the surface charge density increases 
rapidly. It is in this regime that we would expect to 
see most clearly the effects of charge on the dynamics 
of the bubble. Since the maximum radius attained is 
larger, the collapse would be expected to be more violent 
resulting in a higher peak velocity of the bubble wall 
at the moment of collapse and lower minimum radius 
attained. 

The natural frequency of bubble oscillations of small 
amplitude may be found by assuming that the sound field 
can be introduced through a perturbation of amplitude 
a which is small [29| ■ Then by assuming that the bubble 
oscillates about its equilibrium radius Ro, one can express 
its radius at time t as 



R = R (l + x(t)), 



(3) 



where x(t) is a small quantity of order a. Substituting 
Eqn.(3) in the unforced Raylcigh - Plesset equation and 



III. THE BLAKE THRESHOLD 

The growth of a bubble can be determined by various 
threshold conditions [36|]. One is the well-known Blake 
threshold for mechanical growth of a gas bubble. The 
Blake threshold corresponds to the minimum acoustic 
pressure P s = Psiake exceeding which will result in 
the explosive growth of the bubble, culminating in 
cavitation. Blake threshold calculations are made under 
the assumption that the pressure fields are quasistatic in 
nature, and the surface tension dominates over viscous 
and inertial contributions. 

Another threshold that is important is the transient 
cavitation threshold. This is the minimum acoustic 
pressure required for the forced oscillating bubble to 
collapse violently after its maximal expansion at radial 
velocities at least equalling the speed of sound. 
In this section we calculate and discuss the Blake 
cavitation threshold for the bubble. Further discussion 
and results related to the upper transient cavitation 
threshold will be presented in Section V. 

One other threshold condition is that determining the 
diffusion of a gas into the liquid causing bubble expan- 
sion, which is determined by factors including the natural 
resonance frequency of the bubble, the acoustic forcing 
frequency, the gas saturation coefficient, surface tension, 
specific heat ratio and the ambient radius of the bub- 
ble [36|. We do not discuss this threshold, known as the 
rectified diffusion threshold, in this work. 



The pressure of the liquid on the outer surface of the 
bubble wall pi, may be written down: 



p L (R(t))=Pi(t) 
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where the pressure inside the bubble Piit) comprises of 
pressure of the gas and the vapour pressure P v : 
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The net pressure on the bubble wall from the surrounding 
liquid is [32l ]: 



P = [Po-P v 
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where P ex t = Pa +p{t), p(t) being the ultrasound driving 
pressure. 

For our choice of the polytropic index T = 5/3, the 
change in the bubble radius resulting from quasistatic 
changes in the pressure of the liquid Pl (that is, very 
slow pressure changes of pi, with inertial and viscous ef- 
fects being assumed negligible during bubble expansion 
and contraction) outside the bubble may be determined 
from the equation: 



PL 



Pa — Pi, 
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The charge term completely changes the liquid pres- 
sure profile: in particular for small bubbles, the charge 
term dominates over surface tension, reducing its effect. 

For instance for a 5 micron bubble in water at atmo- 
spheric pressure (a fa 0.0725 N/m, Pq = 101 kPa), while 
the surface tension contribution to the pressure in the 
bubble is roughly 2.8 x 10 4 Pa, the effect of introducing 
a small charge of about 0.415 pC would be to reduce 
this by half. Clearly, this has a significant effect on the 
radial mechanical stability of the bubble and the Blake 
threshold that determines the nature of the bubble's ra- 
dial oscillations. The behavior of eqn.(9) is depicted in 
Figure (|3|) for purpose of illustration. In the absence of 
charge, there exists no equilibrium radius below a critical 
value Pq of the pressure - the bubble radius at this point 
undergoes explosive expansion. The presence of even a 
small amount of charge on the bubble surface produces a 
drastic change of behaviour - there exists no equilibrium 
radius for pressures larger than a critical value Pcmax I 
the pressure region Pq < P < Pcmax is a metastable re- 
gion. Our study in this paper is restricted to physically 
realistic regimes, with applied pressures roughly in the 
range 0.4-1.5 bar. 

To obtain the Blake radius for the charged bubble, we 
adapt the procedure for the uncharged case (see for exam- 
ple Harkin et al [37[) to our situation. We first minimize 






FIG. 3: Illustration of the behavior of eqn.(9): P shown as 
a function of R for (top left) without charge and (top right) 
with charge. The introduction of charge changes the plot 
dramatically. The effect of charge Q on a bubble is shown 
schematically for different values of Q, Q2 > Qi, in the bot- 
tom figure. Note that our study in this paper is confined only 
to applied pressures in the range 0.4 x 10 5 — 1.5 x 10 5 Pa. 



equation ([9]) with respect to R, R > 0. This leads to the 
quartic equation: 
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The Blake radius R cr it is given by the real and positive 
root of this equation. We find that: 
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The liquid pressure PL arit corresponding to this critical 
value of the radius is obtained by substituting eqn.(9) 
back into eqn(4) : 



PL crit = Pv 



R 5 
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Rcrit 87rei? crit 



(14) 



The Blake threshold pressure may be obtained from the 
standard definition uHf : 



PBlake — Pq ~ PL C 



(15) 



A rough estimate of the Blake radius and threshold can 
be made for sub-micron sized bubbles for which the con- 
tributions from the static pressure is negligible in com- 
parison with the charge - corrected terms. In this ap- 
proximation, we find that 
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Substituting this in eqn.(8) and using eqn.(9), we find the 
following approximate expression for the Blake threshold: 
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IV. RAYLEIGH COLLAPSE AND THE 

INFLUENCE OF DRIVING PRESSURE 

AMPLITUDE 

After the bubble attains its maximum radius R m ax, it 
proceeds to the main collapse. Its dynamics during this 
phase is described by the Rayleigh equation: 
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(18) 



In considering cavitation in this limit, one is essentially 
considering collapse of a void, ignoring all terms such as 
viscosity, surface tension, etc. The solution for this is 
found to be 



R{t) = Rr 
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(19) 



where at t = t c , the bubble collapses to a point R = 0. 
Rr is a characteristic radius, and T is the time period 
of oscillation of the bubble. It may be noted that the 
characteristic radius R r in this scaling law is different 
from that reported in [38| where the polytropic constant 
was taken to be unity, corresponding to an isothermal 
process. Here we find an estimate for R r for T = 5/3, 
using a similar energy argument as in [38|. Converting 
the potential energy E pot of the bubble at R ma x to kinetic 
energy at Rq we get 
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Using Eqn.(H9]) in (J20j) we obtain 

r-i\ ^ 



R r 



25T^P [An 
2p 



<L 5 (2D 



For T = 5/3, we get 
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The extremely simplified expression, Eqn.(|22|). is none 
the less useful for making some physically relevant ap- 
proximations. This approximation would hold best dur- 
ing the bubble's collapse to a minimum radius R m in- It 
would also be less inaccurate when the dimensions of the 
collapsing bubble arc very small, that is, when R m in is 
very small. This would tend to match more closely those 
cases where the charge on the bubble is high, so that 
the reduction in values of Rmin is correspondingly more. 
We can see that in the presence of a driving frequency ui 
for the system, R(t) would show a frequency dependence 
R(t) ~ w 2 / 5 in the regime near bubble collapse, at higher 
driving pressure amplitudes, so that we have 
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(23) 



a,\ being a prefactor with appropriate dimensions. 
At higher driving frequencies, the bubble typically has 
larger values for its minimal radius, there not being 
sufficient time for complete collapse to occur before the 
expanding regime sets in. Increasing the charge present 
on the bubble enables it to reach smaller dimensions. 
The minimum radius scaled by the driving frequency, 
through R m i n /uj 2 / 5 is shown in the plot of Figure (@Jt). 
As expected from the discussion above, best agreement 
of Eqn. (|2"3")) . as evidenced through a superposition of all 
the curves for different frequencies, is best seen at higher 
charge values and lower R m in values. 

It is to be expected however, that in reality a stable 
bubble can not carry an indefinite magnitude of charge 
Q. This can also be seen on plotting the minimum radius 
Rmin as a function of charge, where, for a given driving 
pressure, the minimum radius of the bubble for every 
driving frequency converges to one value of the charge. 
This, however, needs to be modified as a further physical 
constraint to the system exists in that bubble contraction 
can not also indiscriminately progress indefinitely. The 
smallest dimensions that the bubble can take, that is, the 
least value of the minimum radius R m in reached during 
the bubble's compressive regime, is bounded by the value 
of the van der Waals hard core radius h. The value of h 
is determined by the gas enclosed within the bubble. For 
example, for argon, h has the value h = i?o/8.86. For 
i?o = 5/zm, this equals a value of 0.564 [im. 
Rmin with this constraint is shown in Figure ^p). 
There is a maximal value for the charge a bubble can 
carry for attainment to R m in = h to be possible. This 
limiting value of the charge we denote by Qh- This is 
shown in Figure ^p) , where the minimum radius of col- 
lapse, R m in, has been plotted as a function of charge Q, 
for a bubble of initial radius Rq = 5fj,m, and pressure am- 
plitude P s = 1.35Pq- The minimal bound of R m in = h 



has been shown by a dotted line. The point of intersec- 
tion of the curve for a particular driving frequency with 
this line gives the value of Qh- 
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FIG. 4: (a) (Above): Rmin scaled by uj 2 ' 5 , plotted as a func- 
tion of charge, Q, for various driving frequencies, (b) (Below): 
Minimum radius, Rmin, plotted as a function of charge, Q, 
for various driving frequencies (v = 20, 24, 28, 32kHz); curves 
for higher frequencies are at the top. The horizontal dotted 
line corresponds to the cut-off Rmin = h, the van der Waals 
hard core radius. Ro — 5/im. (Color Online). 

Correspondingly, there is an upper bound on the maxi- 
mum radial velocity of the bubble V max . The presence of 
charge on the bubble serves to reduce the effective magni- 
tude of the surface tension. With increasing charge, the 
bubble is thus able to reach a smaller radius and a greater 
velocity. V max and R m ax as a function of charge Q are 
shown in Fig. ([5j) for high driving pressure of P s = 1.35P). 
The limits to the curves in the plot are due to the limit in 
the maximal value Qh that Q can take for each frequency. 

Furthermore, as the frequency of the driving ultrasonic 
acoustic wave is increased, the maximum radius attained 
by the bubble, R m ax, reduces. This is understandable by 
recalling that the higher the frequency, the shorter is the 
period of negative pressure shear and the bubble is driven 
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FIG. 5: Rmax on the left, and v max on the right, plotted 
as a function of charge, Q, for various driving frequencies 
(y = 20, 22, 24, 26, 28, 30 kHz), lower frequency curves on the 
top. Ro = 5/im, P a = 1.35Pq. (Color Online). 



to cavitation collapse in a shorter time span resulting in 
shorter expansion time of the bubble. This also causes 
the value of the minimum radius to become larger with 
increasing frequency of the driving ultrasound wave. 
The larger the maximum radius reached, the more vi- 
olent the collapse is, resulting in smaller minimum radius. 

The system is very sensitive to changes in the pressure 
conditions. The amplitude P s of the driving pressure 
determines the physically viable minimal radius attain- 
able by the bubble for a given charge. The maximal, 
bounding value of the charge, Qh, reduces progressively 
with increasing amplitude of pressure P s , until beyond 
a critical pressure P m , it is no longer physically possible 
for the bubble radius to contract to such a small value. 
In Figure © we plot Qh as a function of P s for three 
different frequencies (20, 25 and 30 kHz). We show plots 
for two values of ambient bubble radius Ro = 2^tm and 
i?o = 5/iJTi. Each curve demarcates two regions - the 
space below (or to the left of) the curve corresponds 
to the physically permissible region of R m in > h. The 
region above (or to the right of) each curve corresponds 
to R m in < h, which can not be reached in practice by a 
bubble of that corresponding ambient radius. 

It can be seen that the value of the critical pressure 
P m increases with driving frequency, and reduces with 
ambient radius Rq. At very low values of P s , Qh 
becomes essentially independent of frequency for a 
given ambient radius. We denote the pressure where 
this frequency-independence first sets in (approached 
from above) by Pfi. It can be seen that Pfi occurs 
at a lower value Pfi = 1.12Pq (with corresponding 
Qh ~ 1.3 pC) for the larger bubble (Ro = 5fj,m) than 
for the smaller bubble (Pfi = 1.2Po for Rq = 2/x m). 
Moreover there exists a brief crossover region for the 
Ro = 5/i m bubble where the frequency-dependence of 
Qh reverses, before true frequency-independence of Qh 
occurs at Pfi true = 0.9P (for P s > P fl , Q h for higher 
frequencies are greater than for lower values, while for 
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FIG. 6: Plot of Qh as a function of driving pressure am- 
plitude P a , for bubbles of ambient radius Ro = 2/xm and 
5/im. The area below each curve corresponds to the domain 
where Rmin > h. The region above each curve corresponds 
to Rmin < h and is physically unreachable for the bubble of 
that particular Rq. Ro — 5^m, v = 20, 25, 30kHz; 7?o = 2/xm, 
v = 20, 25, 30kHz. (Color Online). 



Pfi true < Ps < Pfi, this is reversed and Q h for higher 
frequencies are lower than that for lower frequencies for 
the 5/xm bubble). 

To understand the effect of the driving pressure, we 
briefly paraphrase below the arguments given in [381 ] for 
the isothermal case adapting it to our adiabatic system. 

Combining the driving sound field with the static pres- 
sure [29| the total external field P ex t can be expressed as 



Pext = Pq{1 - acosuit). 



(24) 



Substituting this in the Rayleigh Plesset equation un- 
der quasistatic conditions: 
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Ro 8neR 4 
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we obtain the quintic equation 
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The behavior of the equation is completely determined 
by the quantity a — 1. For P ext > 0, we have a a stable, 
single solution for R. For negative, small amplitude P ext 
there are two solutions with that at lower R being the 



stable one. 

These two merge only at a critical value P ex t 



P 



Blake 



PBlake < 0. For this, P qas > P: 



P a is always the 



case and leaves the equation without a solution, where 



gas 



0/P 5 and P CT = 2a JR. 



Liquid pressure becoming negative, opposes the con- 
finement effect of the surface-tension contribution, P a . 
Once the bubble is larger than a critical radius R c , pres- 
sure balance at the bubble wall can not be maintained, 
and explosive growth sets in. The bubble is unstable at 
this stage, and further growth leads to increased instabil- 
ity and still more expansion, and quasistatic conditions 
no longer hold. 

With the time period of oscillation T = 2ir/ui being 
larger than the time scale of the bubble's oscillations, 
oscillations of the external pressure can be considered 
as being quasistatic. For crossing the Blake threshold 
Pext < is required so that a > 1. t = gives P ext = 
(a — 1)Pq to be negative. Thus the behavior shown by 
bubbles for a < 1 will be very different from that for 
a > 1; for values of a less than 1, bubble oscillations 
will tend to be less violently expansive and the effect of 
surface tension dominate the dynamics. 



V. EXPANSION - COMPRESSION RATIO AND 
THE TRANSIENT THRESHOLD 

The surface tension greatly influences bubble dynamics 
and can give rise to very distinct behaviours for bubbles 
of different ambient radii Rq, even when all other condi- 
tions are identical. In smaller bubbles, surface tension is 
a very dominant term. Looking at the R m in vs. P s plots 
(Figure 0) for Rq = 2/im and Ro = 5yum, we at once see 
a striking difference between the two. For the larger, 5 
micron bubble, we see that on lowering the pressure P s , 
at a value corresponding to P s = Pfi = 1.12Prj, the R m in 
curves all converge to a point. Lowering P s further brings 
the bubble to a cross-over regime, until reaching a lower 



pressure P s 



J l tr-u 



0.9Pq, below which pressure, the 



curves largely show frequency - and charge-independent 
behavior. As will be seen in the ensuing paragraphs, 
P s = Pfi is actually the transition pressure Pt r beyond 
which violent bubble collapse occurs. 

For the smaller, 2 micron bubble, we do not see any 
cross-over regime, and the pressure P s = Pfi = 1.2P 
where all curves converge such that for P s < Pfi charge 
or frequency-dependence of the curves is suppressed, is 
actually less than the transition pressure Pt r = 1.3Prj for 
Po = 2^m. 

Introduction of charge on the bubble serves to dra- 
matically move P tr to a lower value, with the effective 
surface tension being reduced due to electrostatic inter- 
action, and its effect being enhanced due to the bubble's 
smaller dimensions. 

One obvious measure of the relative extremal values 
of the bubble dimensions is Rmax/Rmin- Other mea- 
sures used to quantify the bubble's dimensions are the 
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FIG. 7: Rmin vs P 3 for different values of w [y = 
20, 25, 30, 35, 40kHz) and Q. Q = OC, O.lpC, OflpC for (a) 
R = 2/mr and (b) Q = OC", OApC, 0.6pC for R = 5/rni. 
(Color Online). 



expansion ratio E = R ma x/ Ro, and the compression ra- 
tio C = R m i n /Ro. Figure ([5]) shows plots of the relative 
extremal bubble radius measure, R m ax/Rmin as a func- 
tion of the amplitude P s of the driving pressure wave for 
two different bubble radii: Rq = 2fi m and R a = 5/j, m. 
We introduce yet another useful and significant measure 
of the relative extent of bubble expansion to compression, 
which we term the expansion-compression ratio, £ 



(=(E-l)/(l-C) = (R ri 



Ro) / (Ro — Rmin) ■ (27) 



Investigating the dependence of this EC ratio ( on the 
amplitude of applied pressure P s yields some delightful 
results and clearly shows the great utility of this mea- 
sure of bubble expansion/contraction. Figure © is a 



FIG. 8: Rmax/Rmin vs P 3 for different values of ui and Q, 
for a) R = 2/im (Q = OpC, O.lpC, 0.2pC); and b) R = 
5/im (Q = OpC, OApC, O.QpC) [u = 20, 25, 30, 35, 40kHz). 
Curves labelled as 1,2,3,4, & 5 correspond to v = 
20, 25, 30, 35, 40kHz, respectively. Increasing Q shifts curves 
upwards & to the left. (Color Online). 



plot of ( as a function of P s for two different ambient 
radius values, Rq = 2/iin and Rq = 5/im. For bub- 
bles with not too small ambient radius Rq (for exam- 
ple, for Rq = 5/im), at low values of P s , the expansion- 
compression ratio £ becomes independent of charge and 
frequency below some P s = P s /, and £ curves for various 
frequencies and charges all superimpose (Figurc(j9jb))). 
This behaviour is not shown by smaller bubbles (for ex- 
ample, for Rq = 2/im), whose £ curves instead show dis- 
tinct charge and frequency dependence even at very low 
amplitudes of driving pressure (Figurc^a))). 

In general for all bubble sizes, the following generic 
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FIG. 9: Expansion-Contraction ratio £ for (a) Ro — 2/im 
(Q = 0pC,0.1pC,0.2pC) and (b) R = 5^m (Q = 
OpC, OApC, 0.6pC)as a function of P s . Curves labelled as 
1,2,3,4, & 5 correspond to v = 20, 25, 30, 35, 40kHz, respec- 
tively. Increasing Q shifts curves downwards & to the left. 
(Color Online). 



behaviour of the EC ratio £ is shown: at lower pressures, 
till a certain pressure P s = P s i, ( shows only very 
weak dependence on charge and driving frequency. In 
Figure (fTUT a)) the dashed and solid curves are shown as 
representative of different charge and frequency values, 
which arc coincident at pressures below P s i. With 
increasing P s , the EC ratio increases to a peak at a 
critical pressure value Pb, followed by a short, steep 



dip upto a second critical pressure value Pt r - This is 
followed further by a regime of an ever-increasing £ with 
increasing P s (Figure (fTtH a))). Between Pb and P tr 
lies a region of negative slope, which is essentially an 
unstable, transient region. 

Surface tension is overcome at a critical radius cor- 
responding to the pressure Pb after which significant 
bubble expansion occurs and the motion is transient 
until pressure Pt r . For Rq being sufficiently large, 
bubble collapse can not be completed fully during the 
compression part of the cycle of applied pressure and 
bubble motion is then stable. This explanation accounts 
for the presence of two thresholds, enclosinga transient 
regime with stable regions on cither side Q. At the 
lower threshold, it can be seen from Figures © QUI) the 
transition from stable to transient conditions occurs 
very steeply 

The lower transient threshold pressure Pb delineates a 
pressure value above which the bubble expansion occurs 
dramatically, and this in fact equals the Blake threshold 
pressure, P b = PBiake- 



It is important to recall at this point that while the 
Blake threshold is a measure of the onset of rapid bubble 
expansion, it g ives us no information at all about bubble 
implosion [39(. The upper transient threshold pressure, 
Ptr, on the other hand, is the value of P s above which 
violent bubble contraction begins. It can be seen, by 
inspecting the R m in versus P s curve (Figure [7]) that at 
P s = P tr , Rmin = Po/2- This, in fact, identifies the 
transition to a strong collapse regime from weaker oscil- 
lations |3J|. This can be verified by incorporating the 
radial velocity into the EC plot. As can be seen in Fig- 
ure (|10T b)). where a scale representing the magnitude of 
maximum radial velocity has been included, it is only at 
driving pressures above the upper transient threshold P tr 
that velocities rise dramatically to high values, where as 
at lower P s , the radial oscillations occur more slowly. 

In Figure (|Tl"T) (a) and (b) we show the frequency re- 
sponse diagrams for the minimum radius and maximum 
velocity attained by a 5 micron bubble in water at 
the low driving pressure of 0.4Po and carrying charges 
Q = 0,0.4,1.0 and 1.24pC. It is clear that increasing 
the magnitude of charge present on the bubble increases 
the magnitude of the response and advances it to lower 
frequencies. The peaks in the response diagram appear 
much earlier, at lower frequencies, for higher charges and 
their enhanced magnitudes (smaller minimum radius and 
larger maximal velocity) point to more violent collapse. 
We can make a further, very rough estimate of the max- 
imal bubble expansion limits at this transient pressure. 
We first consider the following vastly simplifying assump- 
tions: that at R — > R m i„, R — > c, where c is the speed of 
sound in the liquid. We further assume, for purpose of 
this estimation, that as R — ¥ R m ax, R —^ 0. 
At its local minimum at P s = Pt r , ( — £* r will satisfy 
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FIG. 10: (a). The points on the graph where frequency & 
charge -independence sets in (P s i), the Blake threshold (Pi — 
Paiake) and the upper transient threshold pressure (Ptr) are 
shown in the schematic figure. The dashed and solid curves 
typify different charge and frequency values, which superim- 
pose at pressures below P s i. (b). Plot of the £ curve with 
the key for the maximum radial velocity shown on the right, 
for a given frequency and charge. As can be seen, for pres- 
sures greater than the upper transient threshold pressure Pt r , 
bubble collapse occurs violently, at greater velocities. (Color 
Online) . 



{gp )\p,=P tr = 0' so that we get 

dP s 

Now using 

OJXrairi. 



Ro — R)l 



dP s 



= -c 



.ORr, 



dPs 



(28) 



dPs 



2a Q z 



R 5 



4Q 2 



8TreR° 



min 



2a 

min 



R 



dR n 



dP s 

in Equation 
j-tr 



5(P a + f 

Kq 



Q , Rq 



2a 



SneR^R^ 



R 2 



4Q 2 



SnzRmax 



(29) 



, we get 
/ 50 2a 

[r® 

50 



4g 2 



e>2 
^"inax 



^tRliax 



2a 



AQ 2 Ar)c 

rij min lx rain ° " trt min Ix min 



(30) 



where 4>/Rq denotes the equilibrium pressure of the gas 
in the bubble (Eqn.(6)). 

If we further consider the extremal case of R m in ~ h, 
this becomes 
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FIG. 11: Frequency response curves for minimum radius (a) 
and maximum velocity (b) of the bubble. Ro = 5fim, Q = 
OpC, OApC, l.OpC, 1.24pC. (Color Online). 



Now extremizing equation (|31[) with respect to R ma x 
yields a quartic equation for R max — Rm ax at P s = P tr 



max 8nea max 2a ~ U ' 
In the uncharged case, Q = 0, this simplifies to 

1/4 



R 



tr 
max 



/150 
V 2a 



(32) 



(33) 



which is the value of R ma x at the transition point at 
P s = P tr , made under all the simplifying assumptions 
mentioned above. It will be noted that for mid-sized 
microbubbles, for example for Ro = 5/im (plots for 
which have been shown), the point of transition P tr is 
approximately a constant and is largely independent of 
the driving frequency or the charge. 

A comparison of the estimate of R m ax so obtained from 
the above equation for the uncharged case to the value 
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TABLE I: R^ax obtained from Eqn.(32) & graphically. 



Table 1 


Ro 


^max (eqn.) 


RZax (graph) 


2 /jm 


4.87 /mi 


7.14 /mi 


3 /jm 


7.78 /mi 


10.65 /xm 


4 /im 


10.92 /mi 


12.64 /^m 


5 /im 


14.23 /mi 


14.3 /im 


7 /im 


21.3 /im 


17.1 /im 



obtained numerically is given in Table 1 below. Though 
the estimated values are far from accurate in many cases, 
they do, nonetheless, give a quick and useful estimation of 
Rmax at a driving pressure equalling the upper transient 
threshold of pressure P tr - 

These estimates though rather crude, might provide 
useful measures for avoiding undesirable regimes involv- 
ing violent bubble collapses in medical diagnostics and 
applications. 



VI. CONCLUSIONS 

In this work we have investigated the dynamics of 
a bubble forced by an ultrasound field, and seen how 
charge influences its behaviour. Our calculations are for 
a system where an adiabatic equation of state prevails, 
with r = 5/3. We make several interesting observations. 
Charge serves to reduce the effective surface tension of 
the bubble. This causes a charged bubble to not only ex- 
pand to a larger radius as compared to a neutral bubble 
but also collapse to a smaller minimum radius, the lower 
bound of which is given by the van der Waals hard core 
radius. The charged bubble's collapse is also more vio- 
lent, with radial velocities being reached being greater. 
Charge influences and modifies the liquid pressure pro- 
file. The effects are more marked for bubbles of smaller 
dimensions, where surface tension has a predominant in- 
fluence. Studies of the effect of the amplitude of the 
forcing pressure wave show that introduction of charge 
serves to lower the Blake threshold so that the transi- 
tion to violent collapses and oscillations occur at lower 
levels of pressure, especially for microbubbles of smaller 



dimensions and submicron bubbles. We have obtained 
expressions for the Blake threshold in the presence of 
charge. 

We introduce a measure of the extremal dimensions 
reached by a bubble, which is a quantity ( = (E— 1)/(1 — 
C) where E and C are the expansion and compression 
ratios, respectively. The advantage of plotting £ as a 
function of P s is that it captures the distinct positions 
and behaviors of both the Blake threshold pressure as 
well as the upper transient threshold pressure, between 
which points lies a regime of instability for the bubble. 

The maximum magnitude of charge a bubble can carry 
in the system tends to converge to a single asymptotic 
value, Q m ax, and likewise the minimum radius converges 
to a single value for all frequencies, which is however 
modulated by the radial length scale cut-off provided 
by the van der Waals hard core radius. Exceeding the 
magnitude of charge Qh causes the bubble to contract 
to values of R m in that are less than h which cannot 
be physically reached and hence provide a physical up- 
per bound for the system's charge. We also investigated 
Rayleigh collapse for the bubble, obtaining an expression 
for the characteristic Rayleigh radius R r . Frequency- 
dependence of the minimum radius is also captured. We 
also obtain approximate relations for the maximal radius 
at the critical transition point. Most of the numerical 
results presented in this work are for high pressures, cor- 
responding to regimes of violent bubble collapse. We 
have demonstrated the importance of including charge 
in investigating the expansion and contraction of a bub- 
ble under forcing. We have also found scaling relations 
for extremal radial dimensions. Other results, including 
an investigation of the bifurcation structure as a function 
of charge, scaling relations for the maximal charge, etc. 
are being reported elsewhere |34[. 
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